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Abstract 

\^ ", Let (Xt) be a one-dimensional Ornstein-Uhlenbeck process with initial density func- 

tion / : iR + — > which is a regularly varying function with exponent — (1 + rj), r] S 
(0, 1). We prove the existence of a probability measure v with a Lebesgue density, 
depending on r/, such that for every A £ B(M + ): 



o 

C 
3 



lim P/(X t £A\T x >t) = u{A) 



1 Introduction 



Let f2 = C([0, oo),iR) be the space of real continuous functions, and T the standard Borel 
cx-field on Q. For a probability measure /i on (M,B(M)) we denote by P M the probability 
measure on (Q, T) such that B t (w) = w(t) is a Brownian Motion with initial distribution 
/i. If /i = 5a; is the Dirac mass at x G M we denote P x instead of Ps x . Similarly, if /i has a 
density / we use the notation Pf instead of P M . 

ON ■ Consider a one-dimensional diffusion process (X t ) which in differential form may be written 

O i dX t = dB t - a(X t )dt, X = B , (1) 

\Q . 

where the drift a : 1R — >• 1R is a given function. We denote by £ the infinitesimal operator 
^ ■ of the process That is 

Cf := -9xx/ - a9 x /, (2) 



2 

We also denote by C* the formal adjoint operator of C with respect to the Lebesgue measure. 
In other words 

£j: £•/: ir7,„./ • <)Jnf). (3) 

The hitting time of zero T X is defined as the first time that the process (X t ) reaches zero. 
That is Tq := inf {t > : X t = 0}. A similar notation will be used for the hitting time of 
zero for (B t ): T B . 

A probability measure v is said to be a quasi- stationary distribution (qsd) if 

Vt, VA G B(M) P v (X t G A | T X >t) — u(A). (4) 

1 Departamento de Ingenierfa Matematica, Universidad de Chile, Casilla 170, Correo 3, Santiago 3, Chile. 
2 e-mail: Uadser@dim.uchile.cl 
3 e-mail: jsanmart@dim.uchile.cl 



1 



If the drift is regular, this condition is equivalent to the existence of A G M + such that 

Wt,WA G B(M) P v {X t G A, Tq >t) = e- xt v(A), (5) 

which implies that v — v x , the probability measure concentrated on M + with smooth density 
proportional to the unique solution of the differential problem 

\ v^a(O) = , y/ A (0) = 1. [b) 
That is MA G B{M+) we have v x (A) := (/ °° ip x (x)dx)~ l f A <f X {x)dx. 

We remark that in particular from (5) the absorption time is exponentially distributed when 
the initial distribution is u x , that is P Ux (T x > t) = e~ xt . Usually the set of values of A for 
which (6) has a positive and integrable solution is an interval (0, A] and moreover A coincides 
with the ground state of C*. We shall prove this result in the context of a Ornstein-Uhlenbeck 
process, but this holds in many others situations (see [5], [6] and [9] for example). 

This paper deals with the domain of attraction of the qsd. We say that fi is in the domain 
of attraction of the qsd v if 

lim P^X t G • I T X > t) = !/(•), (7) 

t — >oo 

where the limit is taken in the weak topology. We notice that from (4) v belongs to its own 
domain of attraction. We work with absolutely continuous initial distributions fi, whose 
density / is a regularly varying function (see definitions on section 3). Our main result is 
the following theorem, where we assume that the drift is linear a = ax with a > 0, that is 
(X t ) is an Ornstein-Uhlenbeck process. 

Theorem 

Let / : 1R + — > 1R + be a density function with exponent — (1 + rj) where rj G (0, 1). Then for 
all .1 c Bill! ) 

}im P f (X t EA\T x >t) = v x *(A) 

where A* := at] G (0, a). 



We remark that, in the case that the density function f is smooth, the condition lim u {,\' = 
— (1 + rj) is enough to ensure that / has the desired exponent. 

In the literature there are two works directly related with the problem mentioned above. 
The first one was published by P. Mandl [8] who consider general drift assumptions. He 
proved that the bottom of the spectrum of C* is given by A= sup{A G M : tp x does not 
change sign}. Also he proved that under certain hypothesis on the behaviour of the Fourier 
transform of the initial density function / around the point A, the limit in (7) exists and v 
is the probability measure v x . The measure v x is called the minimal qsd. Because of this, 
we say that Mandl's result only deals with the domain of attraction of the minimal qsd. For 
the Ornstein-Uhlenbeck process it is not hard to see that A = a. 
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The other related result is due to S.Martinez et al [9]. They studied the domain of attraction 
of qsd for a Brownian Motion with constant drift a(x) = a with a > 0. In this setting 
A = a 2 / 2. They prove that the limit in (7) exists if Jirn = P > 0. Moreover, when 

(3 G (0, a] v is the probability measure v\* where A* = a(3 — 4^. When (3 > a then v = v\. 

The present paper is organized as follows. In section 2 we present general and well known 
facts about the Ornstein-Uhlenbeck process which we include for the sake of completeness. 
In section 3 we present the proof of the main result of this paper. 



2 Some facts about the Ornstein-Uhlenbeck process 

From now on, (X t ) denotes an Ornstein-Uhlenbeck process, that is to say, the process that 
solves the stochastic differential equation (1) for a linear drift a(x) := ax with a > constant. 

(y — x)'^ 

As usual P x (B t G dy)/dy := ^=e at denotes the transition density of a Brownian Motion 
starting from x. We define the functions h and g by the formulas 

i p —2at p 1at i 

The transition density of the process (X t ) starting from x can be computed as 

p(t,x,y) := P x (X t G dy)/dy = P e - a t x (B h{t) G dy)/dy (9) 

On the other hand, an Ornstein-Uhlenbeck process satisfies a reflection principle. In words, 
conditioning on , using the strong Markov property and the fact that (X t ) and (— X t ) 
have the same law under P , one obtains for x, t > and A G B(M + ) 

P x (X t G A,Tq >t)= P x (X t G A) - P x (-X t G A). (10) 

Consequently, using formulas (9) and (10), and the reflection principle for a Brownian Mo- 
tion, we obtain for every x, t > 0, 

P X (T* >t)= P e -at x (T Q B > h{t)). (11) 

We also have a formula for q the transition density of the submarkovian process (X t l T x >t ) 
given by 



q(t, x, y) = p(t, x, y) - p(t, x, -y) = J^ye ^^'smft \W)J ' ^ 
where it is assumed that xy > 0. 

We remark that formulas (8) and (11) will help us to rewrite probabilities about the process 
(X t ) in terms of the probabilities about the Brownian motion (B t ). Also, before ending the 
section, we notice that the reflection principle for the process (X t ) essentially holds for any 
diffusion process (Y t ) that starting from zero has the same distribution as (—Y t ). Certainly, 
this is the case when a process (Y t ) solves an stochastic differential equation, of the type (1), 
for which the drift a is an odd function and uniqueness in distribution holds. 
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3 Proof of the main result 



Our main theorem relies on the concept of regularly varying functions (for a complete ref- 
erence, see [1], [3]). A non-negative measurable function / : 1R + — > 1R + it is said to be 
regularly varying with exponent (3 (briefly, / has exponent 0), if for all c > 

hm = C P. 

In order to make a discussion in the same terms of our main result, we deal with a function 
/ with exponent of the form —(1 +rj) with rj G M. Since the function f(u)u 1+v varies slowly, 
we have the following asymptotic for In / ( see [1] Proposition 1.3.6) 

hm^M = -(l+v). (13) 
ln-u 

We start the proof of the main Theorem by proving that the set of measures 
{Pf(X t G • | T X > *)} is tight. The next Lemma is technically important for that pur- 
pose. 



Lemma 3.1 

If f : M+ — > M + is a density function with exponent —(1 + rj) for some rj G (0, 1) then for 
all 7 G (0, 77) 

lnn „HJM(# = H <cc , (14) 

u ^°° J xf(x)dx 77 — 7 

Proof: From the hypothesis assumed on / it follows immediately that x 1 f(x) is integrable 
near 00. Therefore from [3], Theorem 1 on section VIII. 9, we have the following limits exist 

hm — — = rj - 7 > 0, 

and 

lim ru (- = 1 - rj > 0, 

u ^°° foxf(x)dx 

from which the result follows. □ 



Lemma 3.2 

If f : M + — > iR + is a density function with exponent —(1 + 77) /or some rj G (0, 1), then the 
set of probability measures: j P/ (A^ G ■ | T X > t) } is tight. 
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Proof: It is enough to check that lim sup Ef (X t 7 | T X > t) < oo for some 1 > 7 > 0, where 

E f [X t | T > tj - j(0jOOj0jOo) , (15) 

and 

:= f B [ D f{x)yiq{t,x,y)dydx, J(A,B,C,D) := (* f° f(x) q(t,x,y) dydx. 

J A JC JA JC 

It can be checked using (12) that there exists a constant c\ = Ci(a) > such that for every 
t > 1 and x E [0, e at ] 

l < J™y~ f q(t,x,y)dy < c 
JoV y <l(t,x,y)dy 
From (15) and the last inequality we have that 

F f Y~< I v fW J(0.e at ,0,°o) f 1 I /(e"«,oo,0,oo) 1 ^ /(0,e af ,0,l) fi , /(e af ,oo,0,oo) 1 

I ^0 ^ L J ^ J(0,e a *,0,l) \ ±_r /(0,e at ,0,oo) J - Cl J(0,e at ,0,l) \ 1 ^ 7(0,e°*,0,l) J 

Therefore to prove the Lemma is enough to find 1 > 7 > such that 

r I^t f(x)y^ q(t,x,y) dydx 

hmsup — < 00. (16) 

Jo Jo f(x)y y q(t,x,y) dydx 

But since < h(l) < h(t) < ^ for every t > 1 and sinh(z) > z for all z > 0, we see from 
(12) that 

re at fl 4 a 3/2 , e «t ( e -^) 2 /•! j,2 

/ / f{x)y~ f q{t,x,y)dydx>^- xf{x)e^^dx y^e^dy. 
Jo Jo y/ire at Jo Jo 

On the other hand q(t,x,y) < p(t,x,y). Thus, using (9) and the last inequality, it follows 
the existence of a constant C2 = 02(a) > such that 

&ff{x)y<q{t,x,y)dydx ^ ^ &at gS/fr) {j^^^d^ ^ 
Jo"' Jo f( x )v y Q(t, x, y) dydx ~~ ~ J eat xf(x)dx 

Now, if < 7 < 1 there exists a constant C3 > such that for x > e at 

poo 

/ y 1 e' a{e ' atx ~ y)2 dy < c 3 ( e - at x)\ 
Jo 

hence from (17) we can conclude that 
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limsup ^/°7 (x)y7g(M ' y) ^ < c.calimsupe^-)^^ 
Jo lo f(x)y^ q{t,x } y) dydx t-oo f e xf(x)dx 



= C2C3 limsup u 



i- 7 C x^f{x)dx 
Joxf(x)dx 



To finish the proof we notice that Lemma 3.1 ensures that the right hand side in the last 
inequality is finite for any < 7 < r\. This proves assertion (16), and therefore the desired 
result follows. □ 

Lemma 3.3 

Let < b < c and x > 0. Then 

p x {T»>b) (ay/* 

- P X (T B > c) ~ Kb) ■ 

Proof: The first inequality is direct. On the other hand (see [7], page 197) after a linear 
substitution we get 

OO 2 OO x 2 

p x (t b >b)_{ I (c _y/> 



P x {T B >c) Ji e -£ du J 1 — i e -£d u 



The last inequality follows since the function u — > (l + ^) is decreasing. □ 
Lemma 3.4 

If f : M + — > M + is a density function with exponent —(1 + rj) where rj G (0, 1), then for 
every t > and c > 

C {U)/C f(x)Pc4T B > t)dx 

— C H/c /W^(T o s >0^ 1 j 

Proof: Since the function x — > P X (T B > t) is increasing on iR + we obtain 

/ f(x)P c 4T B >t)dx<P M u 1 (T B >t) f(x)dx<P M u 1 (T B >t). 

JO U u JO u 

But the function x — > P X (T B > t) is differentiable at x — 0, hence, there exists a constant 
Ci > 0, which depends on t, such that for u > sufficiently large 

f n{U) ' C f(x)P c 4T B > t)dx < c^. (19) 
Jo u u 
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On the other hand, since / is regularly varying, from (13) if k G (rj, 1) we have for large x 
the inequality f(x) > x r^- In particular, for u > large enough we have that 

poo / r\ K r°° I 

L,J^ T » > ? in) L„ > ^ 

Now, since the function u — > is asymptotically decreasing to 0, there exists a constant 
c 2 > 0, which also depends on t, such that for big u > 

[°° f(x)P c *(T* > t)dx > — . (20) 

Jln(u)/c JK ' " V ° ' ~ U K K J 

From (19) and (20), and the fact that k < 1, it follows the result. □ 



Lemma 3.5 

// / : M + — > iR + is a density function with exponent —(1 + rj) where n e (0, 1), £/ien /or 
every s > 

lim P /( T X > t + S ) _ -(ar,)s (21) 

Proof: Let s > and h be the function defined on (8). Then, using (11), it follows for every 
t > that 



P f (T* >t + s) _ J °° f(x)P e - a(t+s)x (T B > h(t + s))d 



x 



P f (T x > t) J °° f(x)P e - atx (T B > h(t))dx ■ 

Notice that h(t) < h(t + s) < ^ and Km h(t) = Therefore, using the bound obtained in 
Lemma 3.3 and setting u = e at , it can be easily seeing that (21) is equivalent to 



Writing 



'e as ln(u) 

Jo° f{x)P e -^^ {Tq > ^)dx \ o e a *in{u) 



V Jn(u)/ 
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from lemma 3.4 and using the substitution y = e as x in the numerator, to check (22) it is 
sufficient to prove that 

Um Q u) f(e-x)P,(T*>± a )d X = e _ a(1+r?)s _ 



u— >oo 



Q u) f(x)P-(T B >±)dx 



But, / is regularly varying therefore lim ^ = e a ( 1+T ') s for all x > 0. From this fact, it 

X'^OO J\ x ) 

follows the result. □ 



Lemma 3.6 



(a) {A | ipx does not change sign } = (— oo,a\; 

(b) For every A G (0, a], / °° ip\(x)dx < oo. 

Proof: Let tp\ be the solution of (6). Then ip\(u) = e au ip\(u) is the unique solution of the 
equation 

C^ x = -A 

<M0) = , ^(0) = 1. 

Hence, to prove 3.6(a), we just need to concentrate our attention on ip\. Note that ip\ is an 
analytic function. Setting ip\{u) = bkU k one sees that b = 0, bi = 1 and one obtains a 

fc>0 

recursion for bk+2 in terms of bk which shows that VA; > &2fc = and for k > 1 

a k t, 1 / \a- 1 \ . . 

hk+1 = j2kTTW.Uo\~^)' (23) 

In particular ip\ cannot change sign if A < a. Notice also that ip a (u) = u. 

On the other hand, if A G (a, 3a) then from (23) we get that b 2 k+i < for every k > 1. 
Thus, lim ip\{u) = — oo. But ip\(0) = and ^(0) > 0, hence, from the last limit, we see 
that there exists x > such that ip\ > on (0, x ) and ip\(x ) = 0. 

Let A > a. We prove then that ip\ has to change its sign. Letting k G (a,min{3a , A}) we 
just proved the existence of some x > such that ip K > on (0, x ) and ip K (x ) = 0. But, 
simultaneously 




Since A > k, by the Sturn-Liouville's theorem (see [12], page 104), there exists y Q G (0, x ) 
such that ip\(yo) = 0. Hence, to prove that ip\ changes sign, it is sufficient to check that 
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il^x{yo) 7^ 0. If this were not the case it can be easily checked out that ip x k \yo) = for all 
k > 0. Thus, ipx = on [0, oo). This can not occur because ip' x (0) ^ 0. Consequently, for 
A > a, ip'xiyo) which implies that ip\ changes its sign. This proves 3.6(a). 



Now, we prove part 3.6(b). For A G (0, a], we have 



f ^ - sg?^ = 2s{ 1+ g2*Tin( 1 

< — J i + y — - — ex P \~\a- 1 y 

- 2a\ & 2k + l P \ &2i + l, 



2i + l 



(24) 



* ^ 1+ S^ e ^£l>>- (25) 



where we have used the fact 1 — z < e 2 , for z G [0,1]. Let d := sup 

fc>0 

d < oo, (25) yields 

1 e2a _ 1 



MA;) - £ 

i=l 



. Since 



The previous inequality shows that f£° ipx(u)du < oo, which proves 3.6(b). □ 
Proof of the Theorem 

Let t' n — * oo. From lemma 3.2, we know that there exits a subsequence t n — > oo and a 
probability measure /x, defined on B(M + ), such that 

j™ p /(^e.|7^>0 = ^ (26) 

Now, the function P X {T^ > s) is continuous and bounded on x G The strong Markov 
property allows us to show that for every n 

Therefore, taking limit as n — > oo and using lemma 3.5, we deduce that for every s > 

POO 

/ P,(T X > s) At (dx) = e-^>> = e- x * s , (27) 
Jo 

where A* := ai] G (0,a). On the other hand, by lemma 3.6(b), we see that y?A* > is 
integrable on [0, oo). This allows us to define for each s > and y > the function 



poo 

A(s,2/):=/ q(s,x,y)(px*(x)dx. 
Jo 



A simple computation shows 

d 

—q(s,x,y) = C x q(s,x,y), 

which is nothing but the classical Kolmogorov's equation. Now, from the given definitions 
we have 



—A(s,y) = C x q(s,x,y)(p x *(x)dx = q(s, x,y)(C* x ip x *)(x)dx 

= -A*/ q(s,x,y)(p x *(x)dx = -\*A(s,y). 
Jo 

Therefore A(s,y) = A(0,y)e~ x * s = (p\*(y)e~ x * s . Hence, integrating this equality over y leads 

roo roc roc 

/ / q(s,x,y)dy<px*(x)dx = / P X (T Q X > s) <p»(x)(dx) = ce~ x * s , (28) 
Jo Jo Jo 

where c = J °° ip x *(x)dx G (0, oo). 

Now, define on B(M + ) the finite signed-measure p(A) := /jl(A) — ^ J A (p x *(x)dx. It follows 
that p([0, oo)) = 0. From (27) and (28), and then (11), we see that for every s > 

PCX) 

/ P e -as x (T B > h{s))p{dx) = 0. (29) 
Jo 



Using the well known distribution for the running maxima of a Brownian motion one has 
P X {T Q B > t) = -j= e^^dy. This identity, a simple substitution and integrating by parts 
allow us to conclude from (29) that for every s > 

/■OO 1 x 2 

/ p[0,x] —= e 2 °Mdx = 0, 

where g is defined as in (8). Finally, because of the set of density functions j^L^e" 2 ^} 

is a complete family and the range of g is [0, oo), we deduce that p[0, x] = 0, for every x > 0. 
Therefore, for all A e B(1R + ), p(A) = ^ f A (p x *(x)dx = v\*{A). Notice that the limiting 
measure does not depend on the initially chosen sequence, from which the result follows. □ 
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